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Abstract 

We propose a unified approach to g-special functions, which are de- 
generations of basic hypergeometric functions 2<fii{a, b; c; q, x). We obtain 
a list of seven different class of g-special functions: i<pi, i<pi, two different 
types of the g-Bessel functions, the g-Hermite- Weber functions, two dif- 
ferent types of the g-Airy functions. We show that there exist a relation 
between two types of the g-Airy functions. 

1 Introduction 

In study of classical special functions, unified theories help us to understand 
special functions. We can study classical special functions denned by differential 
equations easily by confluence of singular points or the method of separation of 
variables. Since we do not have such a unified theory for g-special functions, 
it is difficult to decide whether two different special functions have relations or 
not. For example, we know three different g-Bessel functions and two different 
g-Airy functions and we know one relation between the first and the second 
g-Bessel function. But it is not easy to determine they are completely different 
or not in general. 

Recently many textbooks on g-special functions have been published such 
as Koekoek-Lesky-Swarttouw jT2] , but they have not shown unified theories of 
g-special functions although they list up a huge list of special functions. In this 
paper, we give a unified theory for g-special functions, which come from degen- 
eration of the basic hypergeometric functions 2<pi(a, 6; c; q, x). In our viewpoint, 
we have essentially two different types of g-Bessel functions. We also see that 
two types of g-Airy equations are essentially equivalent but they are different as 
g-series. A connection formula of g-Airy equations is recently found by T. Morita 
[15] . Our list is not enough to study whole of the Askey scheme. It is a future 
problem to expand our unified theory to include the Askey- Wilson polynomials. 

One of the most famous theory on a unified approach to classical special 
functions is by confluence of singularities [TJ [TT] . We denote irregular singular- 
points of the Poincare rank r — I as (r) in Figure [T] We may consider that 
regular singular points are singularity of the Poincare rank 0. We remark that 
we consider the Poincare rank of the Bessel equation at the infinity is 1/2, 
because we think oF±(c; x) is a true Bessel function (see section [2]). The author 
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Figure 1: The coalescent diagram of classical special functions 



lias learned from Professor Hideyuki Majima that the word "confluent" appeared 
in the second edition of Modern Analysis [22] and did not appear in the first 
edition written without Watson in 1902. 

Another approach is by separation of variables of the Laplacian by orthogo- 
nal coordinates. This method is useful to study the Mathieu functions and the 
spheroidal wave functions [TJ [H] . 

In the study of g-special functions, we do not have such a unified approach 
since it is difficult to consider confluence of singular points or separation of 
variables in the g-analysis. 

The third approach is classify differential equations of the Laplace type [2T] 

d 2 y dy 
(do + b Q x)-^ + (ai + blX )^ + (° 2 + b 2 x)y = 

by means of change of variables x — > px + q and y — » g(x)y. We obtain Kummcr, 
Weber, Bessel and Airy functions from equations of the Laplace type. We review 
the third approach in section [2] 

In section [3l we review different types of g-Bessel functions and g-Airy func- 
tions. In this section we introduce two important tools to study g-difference 
equations. One is a shearing transformation and the second is gauge transfor- 
mations by g-products and theta functions. 

The third approach can be easily modified in g-difference equations. We call 
a g-difference equation of the second order with the linear coefficients 

(do + b x)u(xq 2 ) + (ai + b\X)u{xq) + (a 2 + b 2 x)u(x) — 

the hypergeometric type. In section 31 we classify g-diffcrence equations of the 
hypergeometric type and obtain seven types of g-difference equations: 

In section [SJ we compare our classification of g-difference equations of the 
hypergeometric type to special solutions of the g-Painleve equations. 

We may consider g-difference equations in the matrix form Y(xq) — A{x)Y(x). 
The matrix form has more parameters than a single equation of the higher order. 
These parameters are redundant but we can reduce the number of parameters 
since the matrix form of g-difference equation have more transformations than 
the scalar form of higher order. For a g-difference equation of the hypergeomet- 
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Figure 2: The coalescent diagram of q-special functions 



ric type, the corresponding matrix form is a 2 x 2 system 



Y(xq) = (A + A lX )Y(x). 



In section |5] we classify this type of equations. 

The author expresses his thanks to Professor Lucia Di Vizio for fruitful 
discussions when she was invited to the University of Tokyo in the winter 2011. 
The author also thanks Professor Hidetaka Sakai, Mr. Yoshinori Matsumoto 
and Mr. Takeshi Morita for continuous discussions. This work is supported by 
the Mitsubishi foundation and the JSPS Grant-in- Aid for Scientific Research. 

2 Differential equation of the Laplace type 

In this section, we show confluent hypergeometric series and classification of 
second order differential equations of the Laplace type. Although these results 
are already known, we review classical results [21] in order to fix our notations. 

2.1 Confluent hypergeometric series 

We set confluent hypergeometric series found by Kummer. 



Kummer's second confluent hypergeometric series qFi(— ;c;z) satisfies the dif- 
ferential equation 



The infinity is an irregular singular point of the Poincare rank 1/2. The author 
does not know standard term (degenerate confluent hypergeometric series or 
confluent limit hypergeometric series) for qFi. Tricomi |21) used the notation 
Ey(z) = qF\{— ; u + 1; — z)/T(u+l), Jahnke and Emde [7] used another notation 





K(x) = F 1 {--v + l--x 2 /±). 
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It is known that the Bessel functions are represented by confluent hyperge- 
ometric series in two ways. One is related to iF\: 

M») = 2 4(t + l) lFl(1/2 + 1 + 2V] 2%Z) ' 
and the second is related to qFi: 

M*Vz) = r( ^ + 1) o^i(-;^ + i;-^)- 

It is convenient to take qF±(— ; c; z) as a standard form of the Bessel function at 
least theoretically. But in many applications in mathematical physics, 3 v (z) is 
used since it is obtained as the Cylindrical functions. 

2.2 Classification of equations of the Laplace type 

A second order differential equation of the Laplace type 

d 2 y dy 
(a + b ox)-j^ + (ai + b i x )^ + ( a 2 + ° 2 x)y = (1) 



is solved by special functions explained in subsection 12.11 We remark that by 
the transformations (A,/x and h are constants) 

x = \£ + (j,, y = e hx rj, 

•q(£) also satisfy a differential equations of the Laplace type ([TJ. 
We set 

A(h) = a Q h 2 + a x h + a 2 , B(h) = b h 2 + b\h + b 2 . 

Classification of second order differential equations of the Laplace type is as 
follows (see Capitolo Primo in Tricomi's book|21J. 

Theorem 2.1. The equation §Q is solved by special functions iFi(a;c;z), 
J v {z), \Fx{a;\/2\z) and J 1/3 (z). 

(i) Assume that b ^ and A := b\ — 4b b 2 ^= 0. We set 

_ —bi ± b _ a _ A{h) _ ai&o - a bi 

%o ' mhy 11 a ~mhy c ~ 6g ■ 

Then ([T]) is solved by T) — \F\{a\ c; £), where x = A£ + [i and y = e rj. 

(ii) Assume that bo ^ and A = 0. We set 

Then (HJ) is solved by rj = £, a J 2a {l3\f^), where x = A£ + /j and y = e hx rj. 
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(in) Assume that bo = and a$b\ ^ 0. We set 

b? A'(h) A(h) h 

h = ~r> v = — 7 — , a = ^r- > k = -7r- ■ 

b\ b\ lb\ 2ao 

Then Q is solved by r\ = \F\(a\ 1/2; fc£ 2 ), where x = £ + /i one! y = e^^ry. 
Assume that bo = b\ = tm<i ao&2 7^ 0. VFe set 

, ai 4a a 2 — a? , 2 /&7 

" = _ ^ — ' ^ = — ^ — I ' k= v\ ~- 

zao 4ao02 jy % 

TTien ([T]) is solved by rj = \/£</i/3(fc£ 3 / 2 ) 7 where x = £ + [i and y = e hXr q. 
(v) For other cases, (Q} is solved by elementary functions. 



3 g-Hypergeometric functions 
3.1 Basic notations 

We fix the notation of basic hypergeometric series. The q-Pochhammer symbol 
is given by 

(a;q) = 1, 

(a;q) n = (1 - a)(l - aq) ■ ■ ■ (1 - aq'^ 1 ), (n = 1, 2, 3, ...) 



Moo = n^-a^ x )' 

k=0 
n 

(ai,...,a r ;q) n = JJ(aj;g)„. 

i=i 

The basic hypergeometric series is defined by 

1 l l x (ai,...,a r ;g)„ r f")l 1+r_s « 

r V? s (ai,...,a r ; 61,..., b s ;q, z) = > — — — -1 ) n q\*) z n . 

^ (b 1 ,...,b s ;q) n (q;q) n I 1 

The theta function is defined by 

00 

9 q (x):= ]T q n ^/ 2 x n = (q,-x,-q/x;q) 00 , 

n— — oo 

and 9 q (x) satisfies a q-difference equation x9 q {xq) — 9 q {x). 

The theta function 6 q (x) and the q-Pochhammer symbol (a; g)oo are used to 
transform q-diffcrence equations. The following simple lemma plays an impor- 
tant role in this paper. 

Lemma 3.1. We take a q-difference equation 

a(x)u(xq 2 ) + b(x)u(qx) + c{x)u(x) = 0, 
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(1) We set u(x) — v(x)/(sx;q) 00 . Then v(x) satisfies the q-difference equa- 
tion 

(1 — sqx)a(x)v(xq 2 ) + b(x)v(xq) H ^— ^— v(x) = 0. 

1 — sx 

(2) We set u(x) = 6 q (rx)w(x). Then w(x) satisfies the q-difference equation 

^ ^ w(xq 2 ) H b(x)w(xq) + rx 2 c(x)w(x) = 0. 

rq r 

We show two examples: 

(a) For 

a(x)u{xq 2 ) + b(x)u(xq) + (c + dx)u(x) = 0, 
we set u(x) — v(x)/(~dx/c;q) oc . Then we have 

— (c — dqx)a(x)v(xq 2 ) + b(x)v(xq) + cv(x) = 0. 
c 

(b) We set u(x) = 9 q (ax)w(x) for 

axu(xq 2 ) + b(x)u(xq) + cu{x) = 0. 

Then we have 

w(xq 2 ) + b(x)w(xq) + acxw(x) = 0. 
3.2 Shearing transformations 

In the study of differential equations, shearing transformations are useful to 
study irregular singular points whose Poincare rank is non-integer. 

Shearing transformations are also useful for g-differential equations when a 
slope of the Newton diagram is non-integer. 

Definition 3.1. For a q-difference equation 

a(x)u(xq 2 ) + b(x)u(xq) + c(x)u(x) = 0, 
a shearing transformation is the following transformation: 

x = t 2 , p= ^/q, v{t) = u(x). 
The shearing transform of the q-difference equation is given by 

a(t 2 )v(tp 2 ) + b(t 2 )v(tp) + c{t 2 )v(t) = 0. 

We will show an example of a shearing transformation. We take a g-linear 
equation 

(ax + b)u(xq 2 ) + cu(xq) + du(x) = 0. (2) 

Here a ^ 0, b, c, d ^ are constants. The Newton diagram of © at x = oo has 
a slope 1/2. The shearing transform of ([2]) is 

(at 2 + b)v(tp 2 ) + cv(tp) + dv(t) = 0. 
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If b ^ 0, we set w(t) = (p 1 y / — 6/a;p) oc v(i). Then w(t) satisfies 

a(t + y/-b/a)w(tp 2 ) + cw(tp) + d(t - p~ l -b/ a)w{t) = 0. 
If b = 0, we set w(t) — 9 p (t)v(i), as shown in subsection 13. II Then w(t) satisfies 

atw(tp 2 ) + cw(tp) + dtw(t) = 0. 
In any case, w(t) has a regular singularity at t — oo if ad =/= 0. 

3.3 g-Bessel functions and g-Airy functions 

It is known that there exist three types of g-Bessel functions and two types of 
g-Airy functions. It is not clear that relations between such functions. It is 
known that jjp(x;q) and Jv{x;q) are essentially equivalent and J„ (x;q) is 
essentially different. We show that two types of g-Airy functions are related by 
a shearing transformation and Ai q (z) is a special case of Jo [x; q). 

1) g-Bessel functions: 

It is known that there exist three types of g-Bessel functions J^\x; q), J„(x; q) 
and J^\x\q). In most of literatures, Jv(x;q) and jl 2 \x;q) are called Jack- 
son's first and second g-Bessel functions, and Jb {x;q) is called Hahn-Exton's 
g-Bessel function. 

The three g-Bessel functions are defined as 

■ h {x;q)= («?)«» \2) 0(pi [ ' q ;<z ' ^ J' 

(9;g)oo 

The three g-Bessel functions satisfy the following g-difference equations: 
u(a; ? ) - (g^ 2 + q-v^Mxq 1 / 2 ) + u(x) = 0, 

Ji 2) : (l + ^ u(xq) - {q v ' 2 + g" !//2 )n( xq 1 ' 2 ) + u(x) = 0, 

4 3) : u(a:g) + [-(g^ 2 + g^" /2 ) + q u/2+1 x 2 }u{xq 1/2 ) + u(x) = 0. 

It is convenient to take the inverse shearing transform of equations above. We 
set p = g 1 / 2 , t = x 2 and v(t) — u{x). Then we have 

: v(tp 2 ) - [f + p- v )v{tp) + (l + -Q «(t) = 0, 

X (2) : (l + ^) «(^ 2 ) " (P" +P- V )v(tp) + v(t) = 0, 
X (3) : i;(tp 2 ) + +p-") +p» +2 t]v{tp) + v(t) = 0. 



J (1) 
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In the following we set 



El 2) (x;p) 
E ( ^(x;p) 



t" 2 <pi (0,0;p 2u+1 ;p,-x/£) , 
fo*i (-,p 2 " +1 ;p,-x P 2 »+ 2 /4) 
t\ Vl (0;p 2 » +1 ;p,xp 2 /4). 



E)p{x;p) and E_l(x;p) are solutions of the equation Jv X \ E^\x;p) and 

E^l(x;p) are solutions of the equation j'} 2 \ E^\x;p) and E^l(xp~ 2u ;p) are 

solutions of the equation We use the notation E„\x;p) since Tricomi 

used E v {x) in his study in the Bessel functions (see subsection 12. ip . 



By Hahn's formula 



4 2 \x;q) = (~;q) 0O -4 1 \x;q), E<?\x;p) = (-^W^fop) (3) 



In this sense Ji, and Jo are equivalent. Fitouhi and Hamza [2] have defined 
another g-Bessel function j a (x, q 2 ), which is essentially equivalent to 4 (x; q). 

2) g-Airy functions 

It is known that there exist two different types of the g-Airy functions. The 
g-Airy function Ai 9 (x) is found in the study of special solutions of the second 
gr-Painleve equation [9]. The limit q — > 1 tends to the Airy function at least 
around the infinity [5J- The second function A g (x), is called the Ramanujan 
function, which is found by Ismail in the study of asymptotic behavior of the 
g-Hermite polynomial [5J. 



It was not known any relation between the g-Airy function Ai g (x) and the 
Ramanujan function A q (x), but recently Morita has found a connection formula 
between Ai g (x) and A q (x) (see ([8]) in subsection I4.2[) . 

The Airy function is a special case of the modified Bessel function: 



As the same as the differential case, the g-Airy function Ai 9 (x) is related to 



Ai 9 (x) 
A 9 (x) 



i<pi(Q;-q;q, -x), 
<pi(-;Q;q,-qx). 




Jl 6 \x;q). 



Lemma 3.2. If q' 



1, we have 



4 3) {x;q) 



oo 



x u Ai q (—qx 2 ). 



(<?; q) 



00 
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The q-Airy function Ai g (x) satisfies a g-difference equation 

u(xq 2 ) + xu(xq) — u(x) = 0, (4) 

and A q (x) satisfies a q-difference equation 

qxu{xq 2 ) — u(xq) + u(x) = 0. (5) 

These two equation connected by shearing transformation as we has seen in 
subsection 13.21 

Lemma 3.3. If u{x) satisfies (j4]) ; 

v{x) = 9 q (-q 2 x)Ai q (l/x) 
satisfies an inverse shearing transform of a modified equation of Q : 

—q 5 x 2 v(xq 2 ) — v(xq) + v(x) = 0, 
which is is solved by A q 2^x 2 q 3 ). 

We can check out the lemma above directly. 

4 Classification 

We classify a q-difference equation of the hypergeometric type: 

(a + b x)u(xq 2 ) + (ai + b 1 x)u(xq) + (a 2 + b 2 x)u(x) = 0. (6) 

The equation above has transformations which keep the hypergeometric type. 
Such transformations ware known by Hahn [3]. 

4.1 A g-analogue of the Riemann scheme 

Before we list up transformations which keep the hypergeometric type, we in- 
troduce Matsumoto's g-analogue of the Riemann scheme [P3] . 

Definition 4.1. For we set two characteristic polynomials ao/j, 2 +ai[i+a2 = 
and 6o + b±X + 62 =0 The first one is a characteristic polynomial around 
x = and the second one is a characteristic polynomial around x = 00. The 
roots of both polynomials are called characteristic exponents fj,x, fj,2, Ai, A2. The 
characteristic exponents is considered as 00 when ao = or 62 =0. We set 
pi = — aa/bo, P2 = —62/02, which are virtual exponents. 

It is easily checked out that P1P2A1A2/X1/U2 = 1, which is a q-analogue of 
Fuchs' relation for Fuchsian differential equations. If all of exponents are not 
zero, ^ is written only by exponents: 

\i\ 2 (x - px)u(xq 2 ) - {(Ai + A 2 )a; - AiA 2 pi(/ii + p. 2 )}u{xq) + (x - Xi\ 2 pip,2Pi)u(x) = 0, 
which is a g-analogue of Papperitz's differential equation |17j . 
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Definition 4.2. For (JB]), we set a q-analogue of the Riemann scheme: 

!0 oo * 1 
Ml Ai pi ; x > . 
M2 A 2 p 2 J 

The Riemann scheme represents a space of solutions. A g-analogue of the 
Riemann scheme has already shown by Hahn [4], but Hahn's Riemann scheme 
is just a table setting out coefficients a,j,bk- Since exponents is essential in 
differential or difference equations, we array exponents in the scheme as the 
same as the original Riemann scheme. 



We study transformations which keep hypergeometric type. We consider 
two difference equations which are transformed by such transformations are 
equivalent. 



Theorem 4.1. The following four transformations on q-difference equations 
keep the hypergeometric type, (l)x-^-cx, (2)u~^-x 1 u, (3)x^-\/x, (4) u — > 
{p2x;q)oo/(x/piq;q)oou(x). By these transformations, the Riemann scheme is 
transformed as follows: 
(1) x — ¥ cx: 




cx 



= $ 



(2) u — > x 7 it for c = q~< ' : 

oo * 
$ { Pi Ai pi 

k P-2 A 2 P2 

(3) x -> 1/x: 





' 


00 


* 


$ < 


Ml 


Ai 


Pi 




>M2 


A 2 


P-2 



(4) u-> (P2x;q)oo/(x/piq;q)ocu: 

{x/piq; q)c 



oo * 

<I> { p,! Ai px 

P2 A 2 p2 



oo * 
Ai Pi/c 
A 2 cp 2 



oo * 
x 7< f> ^ /ii/c cAi pi 
ui/c cA 2 ,o 2 






oo 


* 


Ai 


/'i 




A 2 


f'2 


Pi 





00 


* 


Ai 


Ml 


M2<T 


A 2 


M2 


piq 



-2 
2 



(M2a;;g)c 



g 2 x 



oo * 
{ pi qjp\P2^\ 1/M2<7 

p 2 q/p-ip-2^2 i/piq 



10 



Remark. If we apply the transformation (4) to the basic hypergeometric series 
2fi (a, b; c; q, x), we obtain Heine's transformation 

(abx/c;q)oo fee ab 
2<pi(a,b;c;q,x) = — -, r\ c; ?, —a; 

We consider the transformations not only for individual solutions but also for 
q-diffcrence equations. Halm's formula ([3]) for the q-Bessel function is essentially 
the same as the transformation (4). When some of the exponents are zero or 
the infinity, we take 9 q (sx) instead of (sx; q)oo in (4) as we have seen on Lemma 

ED 

By means of transformations in Theorem 14. 1[ we easily obtain a part of 
q-analogue of Rummer's twenty-four solutions of the hypergeometric equations. 

Proposition 4.2. The q-difference equation 

(c — abqx)u(xq 2 ) — [c + q — (a + b)qx]u(qx) + q(l — x)u(x) = 0. 

for the basic hypergeometric series 2V?i( a i b; c; q, x) have the following eight so- 
lutions: 

(abx/c^q)^ fee ab 

2(fi{a,b;c;q,x), — 2 <P\ ~, r\ c; q, — x 

(x;q)oo \a b c 

f ^ bq q 2 \ 1 (abx / 'c; q) «, f q q q 2 ab 



1 '2fi\—,—;—;q,x, x 1— — r 2Vi\~,T)—;q,—x 

c c c ) (x;q)oo \a b c c 

aq aq cq \ _ Q (q/x; q) x fq c aq q 

x " 2 <fi [a, — ;— ;q, — , x - — . , . 2^1 7, 7] -r\q, - 



c ' b ' ' a&x/ ' (cq/abx)oo \b' b : b x 

-P (b—-—-n—L\ ~-P (q/ x ''q)oo ^ (q c bq ^ q 



x ^2<fi 0, — ;— ;q, —r- , x "7 — 77 — \— 2^1 1 -,-;—; g, 

\ c a a ox/ (cq/abx;q)oc \a a a x 

Here a — q a ,b — q* 3 and c — q"' '. These solutions are q-analogue of a part of 
Rummer's twenty-four solutions 13f and tends to Rummer's original twenty- 
four solutions when q — >■ 1. 

Proof. We can easily obtained all of eight solutions applying the transformations 
in Theorem 14. II to 2fi(d,b;c;q,x). □ 

Remark. A q-analogue of Kummer's twenty-four solutions are considered by 
Hahn [4]. Other sixteen solutions are not represented by 2^1 ■ We need 2^2 or 
31P2 to represent a q-analogue of Kummer's twenty-four solutions. 

4.2 Classification of g-difference equations 

By the transformations in Theorem 14. 1[ we can classify all of q-difference equa- 
tions ^ of the hypergeometric type. This classification gives a coalescent dia- 
gram of q-special functions. 
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Theorem 4.3. A q-difference equation © of the hypergeometric type reduces 
to one of the following equation by transforms in theorem \4-l\ 

1) When 01036163 ^ , Heine's hypergeometric 2<£>i(a, 6; c; g, 2;): 

(c — abqx)u{xq 2 ) — [c + q — (a + b)qx]u(qx) + g(l — x)u(x) = 0. 

2) When 63 = 0, 01030162 7^ 0, i(pi(a;c;q,x): 

(c — aqx)u(xq 2 ) — (c + g — qx)u{qx) + qu{x) = 0. 

Whenbi = 62 = 0, 03-020163 7^ 0, Jackson's Bessel functions E^ (x;q)F 

u{xq 2 ) - (g' y + q-")u{xq) + (1 + x/4)w(V) = 0. 

3-2) When b\ = 63 = 0, 02 • 030162 7^ 0, Hahn-Exton's Bessel functions 
El 3 \x; q): 

u(xq 2 ) + [-{q v + q- v ) + q 2 - u x]u{xq) + u{x) = 0. 

3- 3) When 63 = a\ = 0, 0262 • 0361 ^ 0, q-Hermite- Weber \<fi{a; 0; q, x) 

axu{xq 2 ) + (1 — x)u{xq) — u(x) = 0. 

4- 1 ) When 61 = 02 = 63 = 0, q-Airy \ip\ (0; — g; g, — x) : 

u(xq 2 ) + xu(xq) — u(x) = 0. 

^-^) When ai = 62 = 63 = 0, the Ramanujan function o<^i (— ; 0; g, —tq): 

qxu(xq 2 ) — u(xq) + u{x) = 0. 

We list up the Riemann scheme corresponding to the classification in the 
Theorem 14.31 In the following, A and B are quasi-constants, i.e. A{qx) = 
A{x),B(qx) = B(x). 

(1) For 2^1(0, 6, ; c; g, x), the Riemann scheme is given by 

!0 00 * 

1 a c/abq ; 

q/c b 1 

The solution space is 

$ = A 2 tpi{a, b;c;q,x) + Bx 1 ~~ t 2 <Pi [—■> 

\c c c 

around x = 0. Here g 7 = c. 

(2) For i<pi (a; c; g; x), the Riemann scheme is given by 

!0 00 * 
1 a c/aq ; x 
q/c 00 
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The solution space is 



aq q 2 



$ = Aiipi(a;c;q,x) + B x 1 7 iV?i — ;— ;q,x 

\ c c 

around x = 0. Here q 1 = c. 

(3-1) The Riemann scheme is given by 

!0 oo * 
q v oo ; x 
q~ v -1/4 

The solution space $ is 

$ = Ax" 2 MO, 0; q 1+ »: q, -x/A) + 5^X0, 0; q 1 ^:^, -x/4) 

around x = 0. 

The standard form of Jackson's first g-Bessel function is 



4 (;g)= (g;g)oo ^^°5« 



T 2 

'+1-0 _lL 

4 



Solutions of the (/-linear equation 



i(xg 2 ) - (<f / 2 + r /2 )«N) + ( x + t) = 



(7) 



is u = Ajjp(x;q) + B J^(x;q). The equation {7} is obtained by a shearing 
transformation of (3-1). 

(3-2) The Riemann scheme is given by 

!0 oo * 
p v oo ; 
p- v oo 

The solution space $ is 

$ = Ax" i<pi (Q;q 1+2v ;q,q 2 x)+Bx- v iy 1 (0; g^ 2 "; q, q^x) . 

(3-3) The Riemann scheme is given by 

!0 oo * 
1 a ; 
oo oo 

The solution space $ is 

$ = Aiipo(a; -; q, x) + B 9(-ax/q)2tpo(q/a, 0; ~;q,ax/q 2 ). 
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(4-1) The Riemann scheme is given by 

{0 oo * 1 
1 oo ; x > . 
-1 oo J 

The solution space $ is 

$ = AAi q (x) + Be xilcix Ai q {-x), 

where lqx = log e xj log e q. 

(4-2) The Riemann scheme is given by 

{0 oo * "1 
1 oo ; x > . 
oo oo J 

The solution space <5 is 

$ = Ao<£i(-;0;<?, -qx) + B 9 q (x) 2 <A>(0, 0; -; q, -x/q). 

The Newton diagram of equations in Theorem 14.31 is as follows. The black 
circle means a coefficient which is not zero. The Newton diagram of (3-2) and 
(4-1) are the same, because the q-Airy function is a special case of E^\x;q) 
in case q v — ±i. The Newton diagram explains that the coalescent diagram in 
Figure [TJ If the Newton diagram of a q-diffcrcncc equation is a subset of the the 
Newton diagram of another q- difference equation, we may take a suitable limit. 



bo bi b 2 
a ai a 2 




a ai a 2 

(1) (2) (3-1) 




(3-2) (3-3) (4-1) (4-2) 



Figure 3: The Newton diagram 



Lemma 4.4. In Theorem \4-3\ is equivalent to (4-1) by shearing transfor- 

mation. 

This lemma shows that there exists a relation between the g-Airy function 
and the Ramanujan function. Recently, Morita |15] has shown a connection 
formula between ki q (x) and A q (x) : 

A q2 (-q 3 /x 2 ) = ^ (-9(-x/q)Ai q (xq 2 ) + 9(x / q)Ai q (~xq 2 )) . (8) 

I9i ii'ZJoo 
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5 Hypergeometric solutions of the g-Painleve equa- 
tions 



As the same as the Painleve differential equations have particular solutions 
represented by (confluent) hypergeometric functions, the g-Painleve equations 
also have special solutions written by g-hypergeometric functions. 

Kajiwara, Masuda, Noumi, Ohta, and Yamada has studied g-hypergeometric 
solutions of the g-Painleve equations [9] . The degeneration diagram of q- hypergeometric 
solutions of the q-Painleve equations is as follows: 

q-P q-Pvi q-P v ^ 9 ~p IV -> q-Pa -> q-Pi 

tt/~i wi ( a ; 0; <7> z ) m \ 

hg a v?i -> ivi -> ^ > ft; z ) -> ivi(0; -«;«>*) -> nonc 

(1) -> (2) -> M -> (4-1) -> nonc 

Comparing our list in Theorem 14. 3[ we do not have (3-1) and (4-2). The equa- 
tion (4-2) is related to (4-1) by a shearing transformation. The equation (3-1) 
appears in another form of q-P\i\. 

It is known that there are several types of the g-Painleve equations, whose 
limit q — > 1 tends to the same Painleve differential equation. It is widely known 
that there exist two different types of q-P\u. One is called q-Pm(A^) by Sakai 
[ZD]: 

yy _ z(z - b 2 t) zz _ y(y - ait) 
a 3 a 4 2-63 h a 4 (y - a 3 ) 

Here the affine root system (Ag ) means the type of the initial value space of 
the g-Painleve equation, which is completely classified by Sakai 20 . 

Another equation is shown by Ramani, Grammaticos and Hietarinta |19j : 

ww (w — ais)(w — SCL2) , 
0304 (w — ai)(w — 04) ' 

which is a symmetric specialization of q-Pyi found by Jimbo and Sakai [8]. 

Jimbo and Sakai found q-Pyi as a connection preserving deformation of a 
linear (/-difference equation. Their q-Pyi(A^) is given by 

yy _ (z - b\t){z - b 2 t) zz _ [y — a\t){y — a^t) ( hb 2 _ ^a x a 2 



a 3 a 4 (2 - 63) (I - 64) 6 3 6 4 (y - a 3 )(y - a 4 ) \b3b4 0304 
For q-PviiA^ 1 ), we set 

t = s 2 , q = p 2 , bi = pen (i = 1, 2), 6, = a* (i = 3, 4), 
y(t) = w(s), z(t) = w(s). 
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Then we obtain the equation (0). In this sense, © as can be considered as a 
symmetric specialization of q-Pyi- Since the type of the initial value space of 
(0) is the same as the type of q-Pyi(A^), we may denote © as q-Piu(A^). 
Kajiwara, Ohta and Satsuma|10j has shown that J^p (x; q) is a special solution 
of q-Pm(A^). This fact is quite natural since J^\x; q) can be represented by 
a basic hypergeometric series 2 <fi(p u+1 / 2 , —p v+x ^;p 2v+x ;p, x), where p = ^/q. 

Theorem 5.1. In the degeneration scheme of q- special junctions in Theorem 
\4-3\ six types of q-special functions appear as special solutions of the q-Painleve 
equations. Type (3-1) is appeared as special solutions of a symmetric specializa- 
tion of q-Pyi. Type (4-2) does not appear as special solutions of the q-Painleve 
equations, but a shearing transform of (4-2) is equivalent to the type (4-1), which 
is appeared as special solitons of q-P\\. 

6 Matrix Form 

In 2 x 2 matrix form, a ^-difference equations of the hypergeometric type is 
given by 

Y(xq) = A(x)Y(x) = (A + Aix) Y(x) (10) 

Here 

Y( X ) = (yf{), a = H a -V Al = fa j»y 

\y2{x)J \a 21 a 22 J \b 2 i b 22 J 

Matrix forms have the following transformations: 

(1) Change x — > cx: A + A\x — » A Q + cA\x 

(2) Change Y -> arT (c = q 1 ): A + Axx -> cA + cA x x 

(3) Change x — > 1/x and Y — > f(x)Y, where f(x) is a scalar function 

(4) Change Y PY: A + A x x ->■ P _1 (Ao + cA x )P 

The equation (|10p has eight parameters Aq and Ax, but we may generically 
deduce to three by the transformations (1-4). 

Lemma 6.1. We assume that b\ 2 = in (|10p . Then y\(x) satisfies a single 
equation of the hypergeometric type. 

Proof. If we eliminate y 2 , we obtain 

y\{xq 2 ) = [an + a 22 + (b n q + b 22 )x]yi(xq) 
+ [(ai 2 a 2 i - ana 22 q) + {a 12 b 21 - a n b 22 - a 22 bn)x - bnb 22 x 2 ]yi(x) = 0. 

By Lemma [331 this equation reduces to the hypergeometric type. □ 

Since we may assume 612 = by a transformation (4), this assumption is 
not essential. We give a classification theorem for (/-difference equation of 2 x 2 
matrix type: 
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Theorem 6.2. A q-difference equation (jlOl) reduces to one of the following 
equation by transforms (1-4) except for the cases (i) Aq or A\ is a zero matrix, 
(ii) det A(x) ee 0. 

1) det A ^ 0,det A\ ^ 0: 2*Pi(a, b; c; q, x) 



A{x) 



1 (l-a)/c\ / -b/c 
q/c J \(c—b)q/c —aq/c 



1 (l-a)/c\ f-l/c 
q/c + \-q/c 



2) det A ^ 0,det Ai ~ diag(0,^), det A{x) ^ const.: i<fii(a; c; q, x) 
A(x) = 

3-1) dctAo 7^0, Af = 1: 

M \ -(<F 1 \ ( 
- ^0 + ^-1/4 

<?-,§,) det A 7^ 0, A x - diag(0,/i),detA(a:) = const : E^\x]q) 



A{x) 



1 o 

-1 q v +q~ u ) 10 -o 2 -" 



a a \ /0 
g — a q — al \0 —q 



3- 3) A ~ diag(0,A), A 1 ~ diag(0,/i): ^(a; 0; q, x) 

A(x) = 

4- 1) detAo ^ 0, A\ ~ diag(0, /i), det A(x) = const.: the q- Airy function 

det^lo ~ diag(0,A),A 2 = 1: £/ie Ramanuj an function 
A(x) ( l ■ ^° °' 

5Mo = M? = i.' 



' 1 o 














i- 




0/ 



Proof. It is easily checked by direct calculation. We show a single equation 
satisfied by y\. 

1) j/i satisfies 

m(^) - c + g " (a + % % 1 M + = o. 
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And y\{x) = 2 fi(a 7 b;c;q,x)/(abx/c;q) OD is a special solution. A fundamental 
solution is 

Y _ 1 / 2<Pi (a, b; c; x) ±5f x 1_7 2 Pi (a?/ c, &g/c; g 2 /c; x) 
~ (a6x/c;g) oc Vy5§ 2 <fii(aq,b; cq; x)x x 1 ~~' 2 ip 1 (aq/c,b/c;q/c;x) 

2) 2/1 satisfies 

2 (c + q- qx) q(c - ax) 

yi(xq yi(qx) H 5 yi{x) = 0, 

c cr 

And 2/i(x) = npi(a; c; g, x)/(ax/c; q)^ is a special solution. A fundamental 
solution is 



1 / np!(a;c;x) ^x 1 ~ f 2 <Pi(aq/c;q 2 /c;x) 

(ax/c;g)oo ^y^i^i (ag; eg; x)x x 1 ~" / 2 <fii{aq/c;q/c;x) 



3-1) 2/1 satisfies 

yi(V) - («" + g-")^i(a?g) + (1 + a:/4)i/i (s) = 0, 

which is the equation of £^'(x; g). A fundamental solution is 

/ x- 2 ^(0,0;g^ +1 ;-x/4) -±— x -* 2(pi (p, 0; g" 2 ^ 1 ; -x/4) N 

V 4(i-^+i) ^ +1 2^i (0, 0; g 2 " +2 i -x/A) x- 2 ^i(0,0;g- 2 --x/4) ^ 

3-2) t/i satisfies 

2/1 (xg 2 ) + [-(g y + q~ v ) + q 2 -" x\ yi {xq) + Vl {x) = 0, 

which is the equation of El (x; q). A fundamental solution is 

Y= (x v up! (0; g 1+2 - g; g 2 x) x~ v m (0; g 1 " 2 - g; g^x) 
^x" i^i (0; g 1+2 ^; g; g 3 x) g-'ar" i^i (0; g 1 "^; g; g 3 - 2y 

3-3) 2/i satisfies 

2/i (xg 2 ) + g(x - 1)2/1 (xg) - aqy^x) = 0, 

And 2/1 (#) = iV?i(a; 0; g, x)/6 q {— ax) is a special solution. A fundamental solu- 
tion is 

y _ 1 f i^i(a;0;g;x) -^6 q {-ax/q) 2 y a (q/a, 0; -; g; ax/g 2 )^ 



x, 



Og(-ax) (a/g;0;g;xg) 6 q (-ax/q) 2 ip (g 2 /a,0; -; g; ax/g 2 ) / 

4-1) 2/1 satisfies 

2/1 (xg 2 ) +x2/i(xg) - 2/1 (x) = 0, 
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which is the (/-Airy equation. A fundamental solution is 



Y 



i<Pi(0;-q;q,-x) e mlcix 1 ip 1 (0;-q;q,x) 
i<Pi(0;-q;q,-xq) -e 7rt ^ x np 1 (0;-q;q,qx)l ' 



4-2) yi satisfies 

yi(xq 2 ) - yi(xq) + x/qyi(x) = 0. 
And yi(x) — A q (x)/9 q (x/q) is a special solution. A fundamental solution is 

Y = 



^ 6 q (l/q) 



o<Pi{-;0;q;-x) 2<Po (0,0;-;q;-x/q 2 ) 



5) yi satisfies 

yi(x(7 2 ) - xy x {x) = 0. 

And yi (x) — l/9 q -2(x) is a special solution. Since A(xq)A(x) = diag(x,xq), this 
case reduces to a (/-difference equation of the first order. A fundamental solution 
is 

(l/9 q ,(x) 9 q2 (x/q)/8 q (x/qT 
\l/e q 2(qx) 9 q2 (x)/9 q (x) 

□ 



Y 



7 Summary 

There exist seven types of (/-hypergeometric equations: 

2<pi, i<pi, 2<^i(0,0;c), i<^i(0;c), i^i(a;0), i^i(0;-g), <£i(-;0). 

Five of seven (/-hypergeometric functions correspond to particular solutions of 
(/-Painleve equations. Jackson's first (/-Bessel 2fi(0, 0;c) corresponds to par- 
ticular solutions of q-Piu, which is a symmetric specialization of q-Pyj. 

In three types of (/-Bessel functions, Hahn-Exton's (/-Bessel jj) (x;q) might 
be a right (/-Bessel function. In two types of (/-Airy functions, Hamamoto- 
Kajiwara-Witte's (/-Airy function Ai q (x) might be a right (/-Airy function. The 
Ramanujan function oVi( — i 0) is connected to the (/-Airy function by a shearing 
transformation and Morita's connection formula. 
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